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ABSTRACT 

In the inner-most regions of radiation pressure supported accretion disks, the turbulent magnetic pressure may 
greatly exceed that of the gas. If this is the case, it is possible for bulk Alfvenic motions driven by the magnetoro- 
tational instability (MRI) to surpass the electron thermal velocity. Bulk rather than thermal Comptonization may 
then be the dominant radiative process which mediates gravitational energy release. For sufficiently large turbu- 
lent stresses, we show that turbulent Comptonization produces a significant contribution to the far-UV and X-ray 
emission of black hole accretion disks. The existence of this spectral component provides a means of obtaining 
direct observational constraints on the nature of the turbulence itself. We describe how this component may affect 
the spectral energy distributions and variability properties of X-ray binaries and active galactic nuclei. 



Subject headings: accretion, accretion disks - 
non-thermal — turbulence - 

1. Introduction 



black hole physics - 
- X-rays:general 



- MHD — radiation mechanisms: 



The thin accretion disk is the model most commonly em- 
ployed in explaining gravitational energy release in active galac- 
tic nuclei (AGN) and X-ray binaries. In the original formulation 
(Shakura & Sunyaev 1973, Novikov & Thorne 1973 NT73), 
both the mechanical structure of the disk and the thermal spec- 
trum which it radiates were solved together, self consistently. 
More recently, the nature of the anomalous "a-viscosity", re- 
quired for angular momentum transport, has been revealed to 
be MHD turbulence driven by the magnetorotational instability 
(MRI, Balbus & Hawley 1991). Spectral models taking into ac- 
count reprocessing of nonthermal X-ray photons, along with the 
implied reflection features (Lightman & White 1988, Guilbert 
& Rees 1988), have proven to be very useful tools in diagnosing 
the variability and geometrical properties of such flows. Fur- 
ther, the idea that the broad-band X-ray power-law continuum 
results from thermal Comptonization of soft seed disk photons 
by a hot diffuse corona (Haardt & Maraschi 1991) qualitatively 
reproduces the observed X-ray continuum in both AGN and X- 
ray binaries. 

Despite these advances, progress towards our understand- 
ing of the mechanics and radiative processes of accretion flows 
have occurred separately with little quantitative connection be- 
tween the two. We propose that bulk Comptonization of soft 
thermal disk photons by the turbulent motions themselves, in 
other words turbulent Comptonization, may substantially al- 
ter the dynamics and the emitted spectrum concurrently. Thus, 
turbulent Comptonization allows for a direct concrete link be- 
tween the details of the disk dynamics and the observed spec- 
trum. 

Turbulent Comptonization is not an entirely new idea. 
Zel'dovich, Illarionov, and Sunyaev (1972) derived a general 
formulation in which turbulent Comptonization could be ap- 
plied within a cosmological context. Thompson (1994) consid- 
ered Alfvenic reconnection-limited turbulent Comptonization 
while calculating the spectral energy distribution (SED) of a 
7-ray burst fireball model. At the same time, Thompson con- 
jectured that this effect may be responsible for producing the 
Comptonized power-law continuum in AGN and X-ray bina- 
ries. Rather than in a corona, we consider here the possibility 
that turbulent Comptonization occurs above the effective pho- 



tosphere, within the body of the disk itself. 

The organization of this paper is as follows. In section 2 
we outline the basic hydrodynamic principles governing tur- 
bulent Comptonization in accretion disks while defining and 
estimating important spectral parameters such as the turbulent 
wave temperature and y-parameter In section 3, we estimate a 
characteristic timescale for turbulent Comptonization from the 
Kompan'eets equation in terms of fundamental thin disk pa- 
rameters. We show the results of Monte Carlo calculations in 
section 4 in order to demonstrate how turbulent Comptoniza- 
tion alters the emergent disk spectrum. We consider how tur- 
bulent Comptonization may observationally manifest itself in 
AGN and X-ray binaries in section 5 and we summarize our 
conclusions in section 6. 

2. Turbulent Wave Temperature and ^-parameter 

In black hole accretion disks, the turbulent velocity of the 
flow may be larger than the thermal velocity of the electrons in 
the inner- most regions. As a result, electrons moving at the tur- 
bulent velocity will be able to Compton up-scatter soft photons 
to higher energy, potentially deforming the observed spectrum. 
Here, we parameterize the turbulent motions in terms of a wave 
temperature T„, which should be compared to the local elec- 
tron thermal temperature of the disk T,, in order to see whether 
turbulent or thermal Comptonization determines the radiative 
transfer At the same time, we deduce an effective turbulent 
y-parameter, allowing us to estimate the contribution to the re- 
sulting non-thermal spectrum. 

2.1. Initial Estimates 

MRI turbulence is fundamentally magnetic in behavior with 
a characteristic velocity given by the Alfven speed va which is 
less than the local sound speed c,. Close to the hole, where 
the radiation pressure is dominant, is given by the radiation 
sound speed, which is larger than the sound speed of the gas. 
Thus, it is possible for the MRI turbulence to be highly super- 
sonic relative to the gas. 

We propose that supersonic MRI eddies in the inner radiation 
pressure dominated region of thin accretion disks Compton up- 
scatter soft thermal photons arising from viscous dissipation. 
In order for this to take place, the electron bulk velocity which 
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corresponds to the turbulent velocity must be greater than the 
local electron thermal velocity. 

To show this, it is convenient to define a turbulent wave tem- 
perature 7^ 



B^w _ i / 2 



(1) 



MRI turbulence is "strong" on the outer scale such that ~ va. 
That is, the kinetic energy density of the turbulence is roughly 

equal to the magnetic energy density (Hawley, Gammie, & Bal- 
bus 1996). If Pw is the magnetic wave pressure given by 



(2) 



we then have a means of relating the turbulent velocities of elec- 
trons to the total (gasH-radiation) pressure of the disk. Recalling 
the fundamental assumption of thin disks, namely the viscous 
stress is related to the total pressure P by the proportionality 
relation ~ aP, we may write for MRI turbulence 

tr<p = (^pvRV^,-^^^'^ = aP. (3) 

That is, the viscous stress is the sum of both the spatially corre- 
lated Reynolds and Maxwell stress of the turbulence. In every 
numerical simulation of accretion disk turbulence, the Maxwell 
stress {BrB ^/4n) exceeds the Reynolds stresses {pvRVcj,) by at 
least a factor of a few. Also as a result of shear, rotation, and its 
magnetic nature, accretion disk turbulence is highly anisotropic 
such that the components of the Alfven velocity obey 

(vl^'j > (var) > (v^j.) . This implies tr^ t^P„, though they are 
of comparable magnitudes. To quantitatively compare tr^ and 
Pw, we introduce a parameter of proportionality ?/' such that 

^ =Pw = tpTR^ = tp(pVRV<l,-^^^\ . (4) 
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ip is a detailed function of the turbulence itself, as it measures 
the level of anisotropy of the flow. Now, we may write the 
square of the Alfven velocity as 



(vl) - 



2Pw , .P 

~ 2aip — 

P P 



or 

(v') ^ <vl)~2aVc?. (5) 

where we have defined P/p = c^. The turbulent wave temper- 
ature is tied to the sound speed and thus the Alfven speed in 
much the same way as the gas temperature is tied to the electron 
thermal velocity. In terms of thin disk parameters, the turbulent 
wave temperature may be written as 



m2 T)2 



(6) 



Here, m is the accretion rate in units of the Eddington accretion 
rate Medd, e is the radiative efficiency, r is the cylindrical radius 
in units of gravitational radii Vg = GM/c^, and A-V me the 
general relativistic correction factors (Riffert & Herold 1995). 

We emphasize that the wave temperature defined above de- 
scribes only the outer scale of the turbulence. Turbulent flows 
may possess power for many decades of wavelength, as is thought 
to be the case with MRI turbulence (Hawley, Gammie, & Bal- 
bus 1996). The general picture of MRI turbulence is that the 



MRI continuously injects energy into the flow at the outer scale 
Ao which is somewhat smaller than the local disk scale height 
H. SmaUer scale modes then remove energy from the MRI near 
the outer scale via mode couplings where the eddy turnover 
time is roughly the dynamical time on the outer scale given 
by the Keplerian orbital frequency When considering a 
scale A smaller than the MRI stirring scale, mode amplitudes 
will most likely be reduced - a point which we shall consider 
later. 

The turbulent wave temperature is independent of black hole 
mass and approaches the electron rest mass value mgC^ jk^ — 
6 X IO'K near the hole for large values of m, a, and for 1. 
Clearly, we see that the RMS bulk motion of the eddies can 
be much larger than the RMS electron thermal motion. For a 
black hole mass of M ~ lO' and V^^Mq, the gas temperature at 
the irmer-edge of the disk is ~ 10' and 10^ K, respectively. 
The turbulent wave temperature, which is independent of black 
hole mass, can be as high as T„ ^ 10** - 10^ K. This is not sur- 
prising since the turbulent stresses, which scale like aP, may 
be up to 10^-10' times greater than the gas pressure at the 
inner edge of black hole accretion disks. Clearly, while mod- 
eUng a thin a-disk the effects of bulk electron motion due to 
vigorous MRI turbulence cannot be ignored in terms of spec- 
tral formation. That is, the contribution to the kinetic energy 
of the electrons originating from random turbulent motions can 
be much larger than the contribution given by random thermal 
motions from which a roughly thermal black body spectrum is 
presumably formed . Notice that T„ cx while the thermal 
black body temperature ~ Jth oc f^l''. Thus, not only will 
be greatest in magnitude close to the inner edge, but the ratio 
7vv/7th will also be largest near the hole. Therefore, effects due 
to large turbulent velocities and wave temperatures will oper- 
ate most dramatically deep within the gravitational potential of 
the hole, where a large fraction of the gravitational power is 
released. Also, disks which extend inward closer to the hole 
possess higher values of T^, leading one to believe that disks 
around rapidly spinning Kerr holes are more prone to turbulent 
Comptonization than say, a disk surrounding a Schwarzschild 
hole. 

2.2. Effects Due to Finite Mean Free Path: Optically Thick 
Reduction 

If the local photon mean free path A^ is smaller than an eddy 
with characteristic size A, the photons are effectively trapped 
by the eddies on that scale. Thus, they will be advected back 
and forth with the flow, not being able to "sample" the turbu- 
lent velocities on the scale A. However, if we presume that 
MRI turbulence possesses power down to the scale of the pho- 
ton mean free path A^, then the radiation may interact with the 
turbulent eddies on scales A < Ap, whose velocity amplitudes 
wiU be reduced with respect to those on the outer scale Aq by 
some amount. This situation is similar to the case of acoustic ra- 
diation in stars where it is approximated that sound waves only 
interact with convective eddies of similar dimension (Goldreich 
& Keeley 1977). 

Consider in particular the example of isotropic Kolmogorov 
turbulence. In this case, the power spectrum of the turbulence 
£(A) (X A^/^, which implies 



V^^.(A) ~ V„(A()) ( 

-^0 



(7) 



where n = 1/3. To calculate the turbulent wave temperature T^, 
on the scale of a the photon mean free path, we make note that 
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T„ oc giving us 



7;,(A;,) ~ 7;(Ao) ( ^ 



2/3 



(8) 



where T^wCAo) is the wave temperature of the outer scale of the 
turbulence given by eq. (6). 

As previously mentioned, accretion disk MRI turbulence is 

anisotropic. In fact, MRI turbulence is completely anisotropic 
such that an eddy of a given velocity amphtude va possesses 
a spatial extent {\r,\^,\^ where \r ^A,^ t^A^. In order to 
proceed without a theory of accretion disk turbulence, we as- 
sume that an inertial range for the turbulence exists such that 
the properties of the turbulence possess a scale-free behavior. 
This allows us to parameterize the velocity amplitude by 



Vw(A)=v(Ar,A^,Aj) ~va 



(9) 



with Li being the dimension of the outer scale eddies in the 
J* direction. In the isotropic hmit, Lr=L^= and r = s = t. 
We must note that accretion disks are Thomson optically thick 
and the diffusion approximation for the radiation transfer is rea- 
sonably accurate within the disk. Thus, a photon scatters 
times before exiting the disk. Since the photon field is roughly 
isotropic, we assume a given photon random walks a distance 
Xp in each direction at a given depth. As a result, photons will 
interact with eddies whose velocities are given by 

We will find it convenient to define 

T]ii=LR/L^ and r]4,=L^/L^ 



(11) 



which allows us to write 



Vw(Ap) : 



A. 



(12) 



where riR,0 > 1 for accretion disk MRI turbulence and n = r+s+ 
t. The corresponding wave temperature evaluated at the photon 
mean free path is then given by 

Tw(Xo) f Xp^ 



T„ {Xp) 



(13) 



2.3. Vertical Structure and y— Parameter 



The vertical structure of accretion disks is one of the biggest 
uncertainties in accretion disk theory. In order to accurately 
model turbulent Comptonization in thin disks, the vertical struc- 
ture must be accounted for. If a soft photon is produced at some 
optical depth tq within the disk, it samples eddies of differ- 
ent scales as it climbs out of the disk, since its mean free path 
changes with density and thus with height. We know neither 
the velocity amplitude of the turbulence nor the density profile 
as a function of depth for accretion disks. In order to move 
forward, we must utilize an ad hoc prescription that takes the 
vertical structure into account in a manner suitable for model- 
ing turbulent Comptonization while not over-constraining the 
disk atmosphere in any particular way. 

One way of dealing with this uncertainty is to vertically aver- 
age the turbulent stresses and the radiative transfer. That is, we 
use the vertically averaged A^ while evaluating the wave tem- 
perature at all depths in the disk. By doing so, we relate Xp to 



the total optical depth t since Xp ~ {nesPT^ and t ~H/Xp. Ef- 
fectively, we are modeling the disk to have a vertically constant 
density and turbulent stress. 

We may relate the outer vertical scale of the turbulence to 
the height of the disk H by introducing yet another parameter 

L, = Ti,H (14) 
where JJz < 1. RecalUng that H/Xp ^ r, the wave temperature 
becomes 



Tw{Xp) : 



-2n 



(15) 



Now, we may write the turbulent y-parameter as 

, . 4^^^^ (Ap) 2 
yw [Xp) ^ t"- 



nigC^ 



yw (Xp) ^ 



OTeC^ riR rj^ nl" 



(16) 



where Vr^ vf* ^ 1 and was chosen since r > 1 for thin 
disks. 

Interestingly, eq. (15) tells us that 7^(Ap) T„{Xo) if 

This is only possible near the inner-edge since 
47rer3/2 BC 

J- ~ ; 

3a m AT> 

for thin disks. Close to the hole, where a large fraction of the 
disk's power is released, A^ approaches Aq as long as m and 
a are sufficiently large. When this is the case, y^^ is roughly 
given by the product of the wave temperature on the outer scale 
T„{Xo) and the square of the optical depth 



(18) 



4kBT„iXo) 2 
2 • 



(19) 
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Fig. 1 . — Turbulent wave temperatures evaluated at the scale A,, using the 
relativistic disk model of NT73 while including Riffert & Herold's (1996) cor- 
rections to the vertical hydrostatic balance. The black holes are maximally 
rotating with a=0.998. Curves (a) and (b) correspond to the disk parameters 
(m,a) equal to (1.0, 1.0) and (1.0, 0.1), respectively. The cascade parameter 
= 2.0 and n = 1/3 for all cases with V'rV'^tIz" ~ ^-^ C') ('-) while 
V'R^vi" = 0.4 for curve (a). Curve (c) represents the case where (m, a) equals 
(0. 1 , 0. 1) with the turbulence confined to the upper layers with an optical thick- 
ness of t/10 . 
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Equations (6) and (18) then allow us to write 



1994, Thompson 1994, and Psaltis & Lamb 1997) 



a V 3 



BC 



(20) 



That is, when \p ^ Ao, the turbulent y-parameter is approxi- 
mately independent of mass and radius. 

Figure 1 shows turbulent wave temperatures for maximally 
spinning black holes of varying m and a. The radius at which 
Th, reaches its maximum value coincides with the most lumi- 
nous radius of the disk, a consequence of T„ directly scaling 
with the accretion stress. The curves shown in Figure 1 are 
independent of black hole mass, a result of Cj and t being inde- 
pendent of black hole mass. The curve labeled "c" depicts tur- 
bulent wave temperatures for m and a = 0.1. For this particular 
case, the turbulence was confined to the upper-most 1/10 of the 
disk. Confinement of the turbulence to a thin upper layer may 
be possible since regions of large magnetic pressure will be rel- 
atively buoyant. Identical wave temperatures and j-parameters 
may be attained for a given m by reducing a by a factor / while 
confining the turbulence to a narrow upper layer whose optical 
thickness is r/f, implying that the local turbulent stresses are / 
times larger than the vertically averaged value. Of course, this 
only applies for our constant density disk atmosphere model 
where each optical depth occupies an equal amount of height. 

The turbulent y-parameters corresponding to the wave tem- 
peratures of figure 1 are shown in Figure 2. For our choice of 
n, oc 1 /r for radii substantially far away from the inner edge. 
Yet, close to the inner edge, and in particular where T„ and y„ 
find their largest values, the dependence of on radius is weak 
which agrees with the estimate given by eq. (20). 

3. CoMPTON Cooling OF MRl Turbulence 

If soft seed photons created by viscous dissipation with en- 
ergy hv ~ ksTe are produced within the disk, they have a chance 
to "resonantly" interact with electrons belonging to eddies of 
spatial scale A ~ Ap. If we assume that the radiation field is 
roughly isotropic, its evolution may be given by the Kompan'eets 
equation (1957) generalized for bulk motion (Hu, Scott, & Silk 
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Fig. 2. — Turbulent y-parameters for relativistic accretion disk models. The 
labels (a)-(c) correspond to the curves in Figure (1) with the same labels. 
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An,y + v^- 



(21) 



n,j is the photon occupation number and tc is the dimensionless 
Compton time where dtc = rteajcdt. Terms only quadratic in 
the turbulent velocity are retained since first order contributions 
vanish when a turbulent eddy is averaged over the domain of the 
flow. It is clear that bulk electron motion produces spectral de- 
formations in much the same way as thermal motion and when 
(v^) > 2>hv/me on the scale Ap, photons stand to gain energy. 

Equation (21) allows us to estimate a characteristic timescale 
?c,w for photons to thermally equiUbrate with the turbulent mo- 
tions. 



d_ 

di 



K,,pC 



3c2 



— ^esPC ^ ■ 



(22) 



\, may also be thought of as the rate in which the eddies cool 
via Compton scattering. In terms of physical disk quantities, 



1 _ W / 2\ 



c Hp c 



(23) 



Where is the Keplerian orbital frequency. Further simpUfica- 
tion is achieved by reaUzing that for thin disks 

^-Oil), (24) 
c 

which is merely a statement thai Ihc disk removes thermal en- 
ergy in the vertical direction by means of radiative diffusion at 
a rate given by ^ (aO)"^ . We are left with 

?-J,~^f22//~r2. (25) 

Note that in the expression above, only eddies on the outer scale 
were addressed and the cascade parameters and 77, were dis- 
carded. As mentioned in section 2.3, A^ approaches Aq near the 
inner-edge for sufficiently large values of rh and a. When this 
is the case for unsaturated turbulent Comptonization, the pho- 
ton spectral index F and the turbulent wave temperature cut-off 
will vary with the turbulence on the dynamical time ' . If the 
turbulent Comptonization is saturated, the photons will be able 
stay in thermal equilibrium with the turbulence, adjusting their 
own "temperature" to match that of the turbulence on the local 
dynamical time of the disk. 

Turbulent Comptonization removes energy from the cascade 
most effectively on the scale \p, while directly placing that 
same amount of energy into the radiation field. The mechan- 
ical energy flux which proceeds down to smaller wavelengths 
is diminished for scales < A^. 

4. Monte Carlo Calculation 

To estimate the effects of turbulent Comptonization on ac- 
cretion disk spectra, we use Monte Carlo radiative transfer cal- 
culations to generate models of integrated emission from the 
innermost radii. The Monte Carlo calculations make use of the 
methods described in Pozdnyakov, SoboF , and Sunyaev (1983). 
We consider a maximally spinning Kerr black hole with emis- 
sion extending to the irmermost stable circular orbit where the 
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stress is assumed to vanish. Since much of the emission then 
comes from material moving at high velocities deep within the 
gravitational potential, we must also account for the effects of 
special and general relativity. To properly model the photon 
geodesies, we use a relativistic transfer code (Agol 1997). 

We calculate the specific intensity emerging from ~ 15 loga- 
rithmically spaced annuli within ~ 20 rg of the black hole. For 
each radius we model the turbulent region of the disk as a ho- 
mogeneous plane-parallel slab. Photons drawn from a Planck 
distribution at the local disk effective temperature Tejj are in- 
jected at the base of the slab and Compton scatter until they 
reach the surface. There is no emission or absorption within 
the slab. For convenience, the scattering electron population is 
drawn from a Maxwell-Boltzmann distribution with tempera- 
ture T„(Xp). MRI turbulence is anisotropic such that electrons 
in an accretion disk will not be well represented by this dis- 
tribution. However if the radiation field is nearly isotropic, 
the Kompan'eets equation (21) indicates that the effects of tur- 
bulent Comptonization depend only on the magnitude of the 
electron velocity, not its direction. The resulting spectra are 
not very sensitive to this assumption since the disk is optically 
thick, implying a highly isotropic photon field. Similarly, any 
velocity distribution with the same second moment will likely 
produce a similar photon spectrum, so our choice of a Maxwell- 
Boltzmann distribution is not restrictive. Further, the luminos- 
ity of a given annulus is normaUzed to the locally available 
gravitational power. To calculate the total spectrum we inte- 
grate over the first 20 Rg using the relativistic transfer code to 
interpolate on our radial grid. For radii where T^jj > T^iXp) we 
assume blackbody emission at T^ff. 

Our Monte Carlo calculations are limited by the uncertainties 
that plague vertical structure models of accretion disks. We ne- 
glect distributed emission and absorption and choose a Planck 
distribution with Tgff for our soft input since it is difficult to 
determine how turbulent Comptonization will alter the verti- 
cal temperature and density profile of the disk. In all cases we 
have calculated, the turbulent slab is effectively thin at the pho- 



ton frequencies of interest, at least at radii at which turbulent 
Comptonization affects the spectrum, so our neglect of absorp- 
tion is justified. Though our model of the disk atmosphere is 
crude, we believe it gives a reasonable estimate of the magni- 
tude and shape of the resulting spectral deformation. 

Figures 3 and 4 show computed spectra for a IQ^Mq black 
hole with different radial profiles of wave temperature and y- 
parameter viewed from an incUnation of 25°. Figure 3 displays 
unsaturated turbulent Comptonization since < 1 for the en- 
tire disk while Figure 4 shows saturated turbulent Comptoniza- 
tion due to this particular disk model's large values of y^. Since 
yw is roughly uniform where most of the X-rays are formed via 
turbulent Comptonization i.e., near the inner-edge, the total X- 
ray spectrum has a well defined slope and spectral index. 

Standard thin disk models predict an effective optical depth 
< 1 for large values of m and a near the inner edge. As a result, 
the gas must reach very high temperatures in order to produce 
sufficient photon flux to balance the turbulent dissipation. Tur- 
bulent Comptonization might alter this picture, because it will 
itself amplify the available soft photon flux as well as prevent 
energy from flowing down to the microscopic dissipation scale. 
The gas temperature may therefore be cooler than would be ex- 
pected in models that only include bremsstrahlung and thermal 
Comptonization. 

5. Observational Implications 

If large turbulent stresses exist in radiation pressure supported 
accretion disks, then bulk Comptonization from the turbulent 
eddies may provide a large fraction of the far-UV and X-ray 
power in black hole or even neutron star sources. The exact 
role turbulent Comptonization will play within a detailed spec- 
tral model is a difficult question to address. The amount of 
power which escapes down to the dissipation scale Xd of the 
flow depends on the details of the turbulence on larger scales. 
At the same time, the properties of the flow on relatively large 
scales are presumably determined by the vertical structure of 
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Fig. 3. — Accretion disk spectrum for a IO^Mq black hole. The turbulent 
Comptonization parameters where chosen from the curves labeled "c" in figs. 
(1) and (2). The corresponding multi-temperature disk blackbody is given by 
the dashed curve. 
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Fig. 4. — Accretion disk spectrum for a IO'^Mq black hole. The turbulent 
Comptonization parameters where chosen from the curves labeled "b" in figs. 
(1) and (2). The corresponding multi-temperature disk black body is given by 
the dashed curve. 
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the disk, which in turn is affected by the amount of mechanical 
energy being converted into thermal energy at A/j. 

In the last decade, paradigms for the broad-band continuum 
spectra for both AGN and X-ray binaries have begun to solidify. 
Most spectral models of black hole accretion disks assume a hot 
diffuse (r < 1) population of thermal electrons above the disk 
(in other words, a corona) in order to reproduce the apparent 
Comptonized spectrum observed in the X-ray band. Two-phase 
models of this sort (e.g. Haardt & Maraschi 1991, 1993) re- 
produce the hard power law tail with the correct spectral index 
r ~ 1 .7 with a Compton reflection hump as long as nearly all 
the gravitational power is locally released in the corona, per- 
haps due to buoyancy and magnetic reconnections of MRI tur- 
bulence (Haardt, Maraschi, & Ghisellini 1994). 

However, none of the current models take into account the 
details of the accretion flow in a predictive manner. If turbulent 
Comptonization substantially contributes to the total output of 
the accretion flow, it then becomes possible to perform both 
spectral and hydrodynamic modeling concurrently since both 
the hydrodynamic and radiative processes are co-dependent. In 
other words, the accretion disk spectrum directly measures both 
the flow, thermal structure, and geometry simultaneously. 

We now discuss how turbulent Comptonization may help ex- 
plain some observations in both galactic and massive black hole 
sources. Our purpose is not to offer definite conclusions, but to 
merely show how this effect may change (or complicate) the 
interpretation of accretion disk spectra. 

5.1. Active Galactic Nuclei 

In many Seyfert 1 galaxies, a powerful extreme ultraviolet/soft 
X-ray feature exists above the extrapolated hard (2-10 keV) 
power law. This component, known as the soft-excess (Wal- 
ter & Fink 1993), is commonly thought to arise from Comp- 
tonization with yet another thermal population of electrons re- 
sulting in a power law that is steeper than that of the 2-10 keV 
continuum. Current disk-corona models cannot account for the 
soft excess in a robust fashion. If the soft excess results from 
Comptonization, it is difficult to explain how and where a "sec- 
ond" corona comes about. A slim disk model (Abramowicz et 
al. 1988) has the ability to produce thermal soft X-rays, but 
the inclination angle of the disk has to be restricted to small 
values. An alternative way of explaining the soft-excess is by 
assuming that it originates from turbulent Comptonization. For 
large accretion rates m ~ 1 and turbulent stresses a ~ 0.1, tur- 
bulent wave temperatures approach T„ ~ 10' K. Lowering rh 
and/or a may yield similar wave temperatures if the turbulent 
stresses are contained within a relatively small number of op- 
tical depths. Narrow Line Seyfert 1 galaxies (NLSl) tend to 
have significant soft excesses which are rapidly varying (Boiler, 
Brandt, & Fink 1996). This impUes that the production of soft 
X-rays takes place close to the hole for these presumably high 
Eddington sources. Turbulent Comptonization may very well 
be the mechanism which produces the soft excesses in NLSl's 
since larger wave temperatures are possible for disks with high 
values of rh while the majority of the power originates close to 
the hole where variability is expected to be most rapid. For ex- 
ample, Figure 3 depicts turbulent Comptonization for a 1O^M0 
black hole accreting at a rate m = 0.1, typical of NLSl sources. 
The spectral index T is greater than two and the thermal cut-off 
is - IkeV. 

Studying the variability properties of NGC 3516, Edelson et 
al. (2000) sampled the optical, UV, soft, and hard X-ray bands 



continuously for nearly three days. They found that the modula- 
tions of the non-thermal X-ray flux and optical/UV thermal flux 
where not correlated. Either up-scattering or reflection mod- 
els require soft thermal photons originating from the optically 
thick disk to either lead or follow variations in the hard opti- 
cally thin component, which apparently is not the case for NGC 
3516. However, this may be reconciled if a special, source- 
dependent, geometry is invoked (Zdziarski, Lubinski, & Smith 
1999 and Chiang & Blaes 2001 in particular, for the case of 
NGC 3516) or if the corona consists of separate discrete blobs 
moving at mildly relativistic speeds away from the disk (Be- 
loborodov 1999). 

A different way of explaining the broad-band variability of 
NGC 35 16 is by assuming the X-rays are produced in the inner- 
portions of the disk by turbulent Comptonization. As long as 
the disk is thin, the hard output will not substantially interact 
with the cooler outer parts of the disk thus temporally discon- 
necting the hard from the soft output. If the turbulent y- param- 
eter is large, then very little UV/optical output emanates from 
the region of the disk that yields hard photons due to the large 
Compton ampUfication. On the other hand, if y^. < 1, then UV 
output which corresponds to the inner X-ray region will vary 
on the thermal timescale fth ^ (a^Y)~^ . From eq. (25), we see 
that the X-ray power in a given band can undergo fluctuations 
on timescale of order O"', which can be much shorter than ^th 
depending on the magnitude of a and the details of the MRI 
cascade. Another possibility is that the soft output from the tur- 
bulent Comptonizing region may contribute relatively little to 
the optical/UV luminosity compared to regions of the disk of 
larger radius as a result of the smaller emitting area. 

In the case of NGC 7469 (Nandra et al. 1998), the lack 
of correlation between the UV and 2-10 keV X-ray luminos- 
ity suggests that neither an up-scattering model nor a reflec- 
tion model can explain the time lags and correlations (or lack 
thereof) between the X-ray and UV output. More recently how- 
ever, Nandra et al. (2000) modeled the RXTE spectral data of 
NGC 7469 and found that the UV power of the source was cor- 
related with the inferred spectral index F, rather than the ampli- 
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Fig. 5. — Flux in the neutral FeK edge band (7.0-7.5 keV) taken from the 
disk spectra computed in fig. (4). The dashed curve corresponds to a disk flux 
(X r~^. 
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tude, of the 2-10 keV band. This implies that the overall X-ray 
spectral feature became softer while the total X-ray luminosity 
increased. Their interpretation of the UV-F correlation was that 
an increase in the UV flux led to enhanced cooling of the hot X- 
ray source which is a direct outcome of Compton up-scattering 
models (see also Chiang 2002). Alternatively, the steepening of 
the X-ray component, along with its increase in power, may re- 
sult from turbulent Comptonization. That is, changes in m pro- 
portionally lead to changes in both the thermal UV luminosity 
and T„. For values of T„ ^ 10^-10^ K, F will roughly increase 
with m as will the soft X-ray luminosity. Further, one would 
expect the thermal cutoff of the high energy X-ray power-law 
to remain approximately fixed if turbulent Comptonization is 
responsible for the enhanced soft X-ray power as along as there 
is no direct physical correlation between the the parameters of 
the turbulent Comptonizing region (7^, y^^) and the electron 
temperature of the hot corona. Coronal up-scattering models 
however, would predict the thermal cut-off to decrease in en- 
ergy. Future missions which have enough sensitivity to detect 
changes in the cut-off energy may therefore be able to distin- 
guish between these models. 

The Fe Ka line emissivity of MGC-6-30-15 may be another 
observational signature of turbulent Comptonization. A signif- 
icant red tail in the fine profile indicates a Kerr geometry (Iwa- 
sawa et al. 1996) where the effects of turbulent Comptoniza- 
tion are potentially at their greatest. Wilms et al. (2001) fit 
the Fe Ka line emissivity to a radially-dependent power law 
oc r~l^ with /3 ~ 4-5. To account for such a steep profile, they 
required a powerful centrally located hard X-ray source with 
incident flux irradiating only the very innermost regions of the 
disk. Those authors as well as others (Ballantyne, Vaughan, 
& Fabian 2003) speculate that the steep emissivity profile may 
result from the extraction of spin energy via magnetic torques 
from the black hole itself (Blandford & Znajek 1977). Perhaps 
a milder way of interpreting the steep X-ray emissivity is by 
invoking turbulent Comptonization. Figure 5 shows the radial 
dependence of the X-ray flux in the 7-7.5 keV band due to 
turbulent Comptonization within a relativistic thin disk. This is 
consistent with a power-law dependency of /? ~ 5. The eddies 
that have the ability to ionize the K- shell of iron are solely lo- 
cated close to the hole, near the marginally stable orbit. This 
provides a natural explanation for concentrating energy release 
in the form of X-rays near the central object. If the Fe Ka line 
emission originates from within the optically thick disk via tur- 
bulent Comptonization, then the ionization and thermal struc- 
ture of the disk must be carefully modeled in order to determine 
whether or not this effect reproduces the observed equivalent 
widths. 

5.2. Galactic Black Holes 

Cyg X- 1 is the most highly studied galactic black hole (GBH) 
source. Like other GBH, it is observed in both the hard (low) 
and soft (high) spectral states (Done 2002) while mostly spend- 
ing its time in the low state. For GBH in general, the soft 
spectral state is dominated by a cool optically thick black body 
while a steep Comptonized power law extends deep into the 
hard X-ray/soft 7-ray region of the spectrum. In the hard state, 
the photon output is dominated by a somewhat flat hard Comp- 
tonized power law with a moderate thermal contribution from 
the cool optically thick disk. The measured luminosities of 
GBH in the hard state are typically lower than the soft state, 
implying a smaUer rh and/or smaUer radiative efficiency com- 



pared to sources in the soft state. 

Gierlihski et al. (1997, hereafter G97) found that in the low 
state, the SED of Cyg X- 1 could not be modeled with a standard 
disk-corona geometry. In order to fit the observed spectra, they 
required a hot 50 keV optically thick (t ^ 6) Wien-like spectral 
component in addition to cool thermal disk emission and a hard 
relatively optically thin (r ~ 1-2) 100 keV power law, pre- 
sumably from thermal Comptonization. They concluded that 
the observed photon- starved spectrum may be explained by a 
geometry similar to the one proposed by Shapiro, Lightman, 
& Eardley (1976) i.e., a geometrically thick hot central corona 
surrounded by a cool optically thick disk. 

In a later study, Gierlihski et al. (1999, hereafter G99) mod- 
eled Cyg X- 1 in the soft state during a hard to soft state transi- 
tion (Zhang et al. 1997). The hard power law extended to ~ 600 
keV with a spectral index F ~ 2.5 implying a very optically thin 
plasma with r ~ 10"^. If thermal Comptonization were respon- 
sible for the hard continuum, bumps would appear in the X-ray 
spectrum due to the distinct scattering orders. However, these 
features are not present. As a solution, G99 proposed that the 
Comptonizing corona was composed of a hybrid thermal/non 
thermal plasma whose distribution function was roughly the 
sum of a relatively cool ~ 50 keV MaxweUian and a power 
law. 

Zhang et al. (1997) noted that during the state transition, the 
bolometric luminosity of Cyg X-1 changed very Uttle. From 
fliis fact, Poutanen, Krohk, and Ryde (1997) PKR97 put forth 
the proposition that the observed state transition of Cyg X-1 
corresponded to a change in the mechanical state of the disk 
characterized by the inner edge moving inward as the spectral 
state changed from hard to soft. This claim was supported by 
the fact that the temperature of the thermal disk component in- 
creased as Cyg X-1 evolved into the soft state. The increase 
in the number of soft photons then has the ability to cool the 
corona, thus steepening the hard continuum. 

Adopting this general picture, further details regarding the 
accretion flow of Cyg X-1 can be inferred if the ~ 50keV spec- 
tral component fitted by G97 and G99 is assumed to result from 
turbulent Comptonization. If we approximate that the bolomet- 
ric luminosity L = eMc^ and are constants during the state 
transition, then >',,, is the relevant parameter which is allowed to 
change. In the hard state, the 50 keV component was fit with a 
Wien-law, implying saturated Comptonization while in the soft 
state, the 50 keV component was fit with a steep power law, typ- 
ical of unsaturated Comptonization. Thus, during the transition 
from the hard to soft spectral state, y^, decreased. As previously 
discussed, the outer scale of the turbulence Aq may approach 
the photon mean free path Ap near the inner edge of the disk, 
where the spectral contribution from turbulent Comptonization 
is greatest. With this, eq. (20) shows us that y„ oc a"'. Thus, 
a increased as y^ decreased during the hard to soft transi- 
tion. This is consistent with the constant luminosity and wave 
temperature condition by noting the dependencies of the wave 
temperature on the outer scale given by eq. (6). That is, the 
quantity lii^ /e^r^ must decrease during the transition. A con- 
stant luminosity requires rh oc e"' which constrains 1 /e^r^ to 
decrease. But, e oc yielding the result that r must have de- 
creased, which is consistent with the hypothesis of PKR97. 

To summarize, by making the assumption that the 50 keV 
component in both the hard and soft spectral state of Cyg X-1 is 
due to turbulent Comptonization we may then determine further 
details of the accretion flow. If the proposition of PKR97 is 
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utilized to explain the spectral evolution of Cyg X-1, the basic 
model dependencies of turbulent Comptonization necessarily 
require that a increased while m decreased during the spectral 
transition. 



6. Summary and Conclusions 

Turbulent Comptonization may prove to be an important me- 
diator of gravitational energy release in accretion disks. If so, 
simultaneously modeling both the mechanics and radiative pro- 
cesses of the flow not only becomes possible, but necessary in 
order to accurately fit the data. We have shown that the ampli- 
tude and shape of the spectral feature resulting from turbulent 
Comptonization coincides with observed components in both 
AGN and X-ray binaries. If turbulent Comptonization does not 
manifest itself spectrally, then its absence in the SED of ac- 
cretion flows constrains values of fundamental disk parameters 
such as a,m, ande. 

We emphasize that throughout this paper, we have assumed 
that turbulent accretion stresses may greatly exceed the gas pres- 
sure. It is this fact that allows turbulent Comptonization to op- 
erate as an important process. It has been argued that in ra- 
diation pressure supported disks, buoyancy limits the turbulent 
magnetic pressure to values below that of the gas (Sakimoto & 
Coroniti 1981, 1989; Stefla & Rosner 1984). Further, we as- 
sumed that large turbulent stresses exist down to the scale of 
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